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abstract’^ 

The  ten  conditions  which  must  be  satisfied  before  simple 
magnetohydrodynamic  (MHD)  waves  can  propagate  in  a  plasma  are 
enumerated  and  discussed.  Assuming  that  all  but  the  tenth  condi¬ 
tion  are  satisfied  throughout  the  ionosphere  and  that  the  earth  is  a 
perfectly  conducting  plane,  the  propagation  of  MHD  waves  in  a 
horizontally  stratified  ionosphere  composed  of  electrons,  ions, 
and  neutrals  is  studied  when  an  arbitrarily  oriented  plane  wave  is 
obliquely  incident  at  the  lower  edge  of  the  ionosphere.  The  equations 
governing  the  behavior  of  the  field  quantities  are  cast  in  the  "coupled 
form"  suitable  for  computation  by  means  of  a  high-speed  modern 
computer,  but  no  numerical  results  are  obtained  because  of  the 
prohibitive  amount  of  work  necessary  to  solve  the  equations  in  the 
case  of  oblique  incidence. 

The  velocity  of  MHD  waves  is  =  B/’\l  [App.  In  a  fully 
ionized  gas,  p  is  the  mass  density  of  the  plasma:  p  =  Pe  ■*"  Pi* 
is  shown  that  in  a  partially  ionized  gas  p  is  the  effective  density 
of  the  material  responsible  for  the  waves  and  that  for  sufficiently 
low  wave  frequency  and  for  a  weakly  ionized  plasma,  the  mass 
density  of  the  neutrals  should  be  included  in  p  or  p  =  Pe  Pi  Pn* 
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I.  INTRODUCTION 

Recently,  the  propagation  of  magnetohydrodynamic  (MHD)  or 
hydromagnetic  (HM)  waves  in  the  ionosphere  has  been  studied  extensively 
by  physicists  and  geophysicists  alike.  Physicists  are  interested  in  this 
phenomenon  as  an  example  of  the  propagation  of  very  low  frequency 
electromagnetic  waves  in  an  inhomogeneous,  anisotropic,  and  bounded 
medium,  1-0  while  geophysicists  are  interested  in  it  as  a  possible 
explanation  of  geomagnetic  micropulsations,  ^  giant  pulsations,  ®  ionospheric 
noises,  the  heating  of  the  ionosphere,  9-11  and  the  propagation  of  HM 
pulses  that  are  generated  by  high-altitude  nuclear  explosions.  12 

Several  conditions  must  be  satisfied  before  MHD  waves  can 
propagate  in  a  partially  ionized  plasma  permeated  by  an  external  mag¬ 
netic  field.  The  more  important  ones  are*  (a)  the  wave  frequency  gj 
must  be  much  less  than  the  ion  cyclotron  frequency  (b)  the 

phase  velocity  Vpj^  =  u/k  and  the  Alfvefi  speed  must  both  be 
much  less  than  the  velocity  of  light  c;  and  (c)  the  conductivity 
must  be  high  enough  so  that  the  inequality  u/(po<r)<5CV?  is  satisfied.  14 
In  addition,  if  the  following  conditions  are  also  satisfied  the  theory  of 
HM  waves  is  relatively  simple:  (d)  the  static  magnetic  field  is  uniform 
and  satisfies  the  condition  for  appreciable  HM  coupling,  15  i,  e.  , 


(1) 


(e)  the  pressure  and  gravitational  forces  are  negligible  compared  with 
the  electromagnetic  forces;  (f)  the  plasma  is  in  thermal  equilibrium 
and  is  inviscid;  (g)  collisions  of  the  charged  species  with  themselves 
and  with  each  other  are  negligible  compared  with  collisions  between  the 
charged  species  and  the  neutrals  because  the  plasma  is  so  slightly 
ionizedlo.  1';  (h)  the  velocity  of  the  neutrals  can  be  neglected;  (i)  the 
Alfvdh  speed  Va  is  much  greater  than  the  sound  speed  Cq;  and 
(j)  the  electron  collision  frequency  is  much  greater  than  the  electron 
cyclotron  frequency  Ug  and  the  ion  collision  frequency  Vj  is  much 
greater  than  the  ion  cyclotron  frequency  uj. 

A  few  comments  on  the  domain  of  validity  and  the  physical 
significance  of  these  ten  conditions  are  appropriate.  Condition  (a) 
is,  in  fact,  a  definition  of  MHD  waves  as  plasma  waves  whose 


See  the  Glossary  for  the  definition  of  the  more  important  symbols  used 
in  this  report. 
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frequency  is  much  less  than  the  ion  cyclotron  frequency.  When  Wi, 

the  ions  cannot  be  considered  to  be  stationary  as  they  can  in  the  opposite 
limit  when  When  (b)  holds,  the  displacement  current  may  be 

neglected  in  comparison  with  the  conduction  current  in  Ampere's  Law 
V  X  =  J  +  (85/8t)  “  7,  and  the  use  of  nonrelativistic  MHD  equations 
is  correct.  Physically  (b)  means  that  the  rest  energy  density  of  the  plasma 
pc^/2  is  much  greater  than  the  energy  density  of  the  magnetic  field 
B^/(2|1q).  The  e -folding  distance  d  is  the  distance  in  which  the  amplitude 
of  the  MHD  waves  is  attenuated  to  e"^  =  0.  368  of  its  initial  value.  It 

is  easily  shown^^  that  in  a  plasma  of  large  but  finite  conductivity 


d  = 


Zuocrv 
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In  order  to  ensure  that  the  MHD  waves  are  only  slightly  damped, 
condition  (c)  must  hold  sodiat  d/X  =  d/X^^  1,  the 

waves  are  heavily  damped  in  a  fraction  of  a  wavelength,  and  it  makes  no 
sense  to  talk  about  wave  propagation  in  the  plasma.  Likewise,  if  we 
assume  that  the  amplitude  of  the  MHD  waves  varies  with  time  as  exp(jut) 
and  allow  u  to  be  complex,  we  must  require  that  the  time  in  which  the 
amplitude  of  the  waves  is  reduced  by  a  factor  of  l/e  of  their  initial 
value  be  much  larger  than  one  period  of  oscillation  of  the  waves.  It 
can  be  shown  without  much  difficulty  that  the  latter  condition  requires 
that  the  wavelength  of  the  MHD  waves  be  much  greater  than  a  critical 
wavelength  Xc  =  tt/p.oVao'-  IMc)  is  satisfied,  then  X^f§>X^.  Condition 
(d)  guarantees  that  the  coupling  between  the  plasma  and  the  static  mag¬ 
netic  field  is  strong,  i.  e.  .  that  the  dynamic  and  magnetic  forces  are 
comparable.  (If  the  coupling  is  weak,  we  are  in  the  hydr_odynamical 
regime.  )  The  requirement  that  the  seed  magnetic  field  B  be  uniform 
is  not  essential,  but  matters  are  simplified  greatly  when  it  is  true. 

The  magnitude  of  the  earth's  magnetic  field  can  certainly  be  assumed 
to  be  constant  throughout  the  entire  ionosphere;  however,  its  changing 
direction  should  be  taken  into  account.  Like  (a),  we  can  regard  con¬ 
ditions  (e)  and  (f)  as  those  necessary  for  pure  MHD  waves  to  be 
excited  in  a  plasma.  The  inclusion  of  pressure,  gravitational,  and 
dissipative  forces  adds  nothing  important  to  the  picture,  but  their 
omission  results  in  great  simplification.  Conditions  (g)  and  (h)  are 
somewhat  more  difficult  to  interpret  physically;  a  good  discussion  of 
them  is  given  by  Fejer.  Investigators  in  controlled  fusion  call  con¬ 
dition  (i)  the  "low-beta  condition,"  P  being  the  ratio  of  the  sum  of  the 
gas  and  plasma  pressure  to  the  magnetic  pressure.  Thus, 
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PoP 


The  magnetic  pressure  is  approximately  0,01  dyne/cm^  throughout 
most  of  the  ionosphere.  At  130  km  altitude,  the  sum  of  the  neutral  gas 
and  plasma  pressures  amounts  to  about  0.  01  dyne/cm^  and  decreases 
rapidly  with  height  so  that  it  is  only  10-4  dynes/cm^  in  the  Fz  layer, 
or  about  300  km.  It  follows  that  the  assumption  that  condition  (i)  holds 
is  reasonable.  From  (12)  it  can  be  easily  shown  that  when  condition 
(j)  holds,  the  use  of  a  scalar  conductivity  instead  of  a  tensor  conductivity 
is  justified.  In  his  original  paper,  in  which  the  existence  of  MHD  waves 
is  predicted,  and  in  his  book,  Alfve'n  used  a  scalar  conductivity.  Con¬ 
dition  (j)  shows  when  this  is  permissible. 

In  this  report,  we  shall  assume  that  all  of  the  ten  conditions 
hold  except  condition  (j).  The  ionosphere  will  be  assumed  to  be  a 
horizontally  stratified,  three-component  mixture  of  electrons,  ions, 
and  neutrals  whose  properties  vary  with  altitude  as  shown  in  Table  I. 
From  a  study  of  Table  I,  it  can  be  seen  that  the  ionosphere  is  a  weakly 
ionized  plasma  and  that  condition  (g)  holds.  Further,  it  is  assumed 
that  the  earth  is  a  flat,  perfectly  conducting  medium  at  the  low  fre¬ 
quencies  in  which  we  are  interested,  viz,  .  f  <  1  cps.  The  boundary 
between  the  bottom  of  the  ionosphere  and  the  upper  atmosphere  and 
the  top  of  the  ionosphere  and  the  exosphere  will  be  sharp  if  the  wave¬ 
length  of  the  MHD  waves  X.  is  much  larger  than  the  dimensions  over 
which  the  plasma  density  decreases  to  zero.  Following  Spitzer^^  and 
estimating  the  characteristic  length  L  by  twice  the  scale  height,  we 
conclude  that  the  lower  boundary  of  the  ionosphere  may  be  regarded 
as  sharp,  while  the  top  boundary  must  be  considered  to  be  diffuse. 

The  latter  conclusions  follow  from  the  fact  that  LgQ  km  =  km 

and  L550  km  -  200  km,  and  for  f  =  1  cps,  \  ^  100  km,  (For 
estimates  of  the  scale  height  see  reference  20.  ) 

In  a  recent  paper,  Francis  and  Karplus,  using  the  above 
model  and  making  the  above  assumptions,  studied  the  propagation  of 
HM  waves  in  the  ionosphere.  However,  they  confined  themsleves 
to  a  vertically  incident  plane  monochromatic  wave  near  45  degrees 
geomagnetic  lattitude.  Here  we  should  like  to  treat  the  more  general 
case  of  an  obliquely  incident  wave  and  indicate  in  what  sense  the 
problem  can  be  considered  to  be  solved.  No  numerical  results  have 
been  obtained  with  our  equations  describing  oblique  incidence  because  of 
the  large  amount  of  machine  time  necessary  to  solve  them.  It  is 
believed,  however,  that  with  the  help  of  the  approach  detailed  in 
Section  II,  useful  numerical  results  can  be  obtained  if  one  has  un¬ 
limited  financial  resources. 
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TABLE 

■"<“  :=r,i 

N  X  lom  •  ^  numbers  expres^sed  as  Nm  signify 


Altitude, 

km 

T.  Ok 

Particles/ 

cm^ 

Ni  =  Ng 

lons/cm^ 

W 

<^e  ‘'i 

V 

e 

f 

n 

80 

205 

4.  30^4 

1.  oo3 

29.  0 

1.60^  8.47®  2.05® 

3.  29® 

r 

2.8® 

90 

225 

6.  70*^ 

2.  50^ 

29.0 

1.  60^  8.43®  3.  24^ 

5.  42® 

1.4® 

100 

280 

1.1013 

1.  22® 

28.4 

1.  62^  8.  39®  5.  37® 

1.00® 

3.  2® 

120 

435 

6.20II 

1.90® 

26.  7 

1.  71^  8.  32®  3.  12^ 

7.  93® 

3.  9® 

140 

620 

I.45II 

2.355 

25.  7 

1.  76^  8.24®  7.44I 

2.  66^ 

4.3‘ 

160 

780 

5.40IO 

2.  55® 

24.  7 

1.81^  8.  16®  2.  83I 

1.  37® 

4.  S*’ 

180 

905 

2.  50l0 

2.  70® 

23.  7 

1.  87^  8.09®  1.  33I 

8.86^ 

4.  6® 

200 

995 

1.  2810 

2.  75® 

23.  0 

1.91^  8.02®  6.94O 

6.  462 

4.  7® 

250 

1190 

3,  54^ 

3.75® 

21. 0 

2.042  7.84®  2.O1O 

3.  97^ 

5.  5® 

300 

1300 

1. 40^ 

3.  35® 

19.  5 

2.  162  7  ^^6  8.  25-1 

3.832 

5.  2® 

350 

1315 

6.  44® 

3.  90® 

18.  4 

2.  232  7.49®  3.  91*1 

4.  20^ 

5.6^ 

400 

1320 

3.  10® 

4.  22® 

17.  4 

2.  3i2  7.  336  2.  94-1 

4  432 

5.8® 

450 

1325 

1.  58® 

4.  155 

16.  8 

2.  342  7.  17®  1.  00*1 

4.  32^ 

5.  7® 

500 

1325 

8.  10^ 

4.00® 

16.  2 

2.  372  7.01®  5.  26*2  ^ 

4.  14^ 

5.65® 

550 

- 

1325 

4.  377 

3.  60® 

16.  1 

2.  342  6.  86®  2.  86*2  : 

1.  732 

5.  35® 

W.  E.  Francis  and  R.  Karplus,  loc.  cit. 
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According  to  the  simple  theory  of  HM  waves  in  a  fully  ionized 
gas,  the  velocity  of  HM  waves  isl4 


(4) 


In  a  partially  ionized  gas.  the  question  naturally  arises:  "What  is  p 
in  eq.  (4)?"  In  Section  HI.  it  will  be  shown  that  p  is  the  effective 
density  of  the  material  responsible  for  the  HM  waves.  If  the  frequency 
of  the  waves  is  low  enough  and  the  ion  collision  frequency  is  high 
enough,  the  velocity  of  the  ions  may  be  synchronized  with  the  velocity 
of  the  neutrals  as  indicated  by  (38).  _This_occur8  when  Q  =  2wNn/(NiV;) 
is  small  compared  with  one.  Then  v^  ^  v^,  and  we  must  put 
p  =  Pe  +  Pi  +  Pn  (4).  As  a  result,  V;\^in  a  weakly  ionized  plasma  can  be 
much  less  than  its  value  in  a  fully  ionized  plasma.  On  this  basis,  the 
rapid  decrease  of  in  the  D  and  E  layers  of  the  ionosphere,  as 
shown  in  Fig.  1.  is  readily  understandable.  The  magnitude  of  Va 
is  of  considerable  physical  significance  because  the  e-folding  distance 
d  depends  on  it  as  shown  in  (2). 


II.  THE  PROPAGATION  OF  PLANE  MHD  WAVES  OBLIQUELY 
INCIDENT  AT  THE  LOWER  EDGE  OF  THE  IONOSPHERE 

In  this  section,  the  equations  governing  the  propagation  of  MHD 
waves  in  the  ionosphere  will  be  enumerated  and  used  to  determine,  in 
principle  to  be  sure,  the  reflection  and  transmission  coefficients  of  a 
plane  MHD  wave  obliquely  incident  on  a  horizontally  stratified 
ionosphere.  Conditions  (a)  to  (i),  inclusive,  will  be  assumed  to  be 
satisfied  throughout  the  ionosphere.  In  truth,  there  are  regions  of 
the  ionosphere  where  conditions  (c)  and  (d),  as  expressed  by  d/X^l 
and  (I),  are  not  satisfied  and  there  the  quasi -hydrodynamical 
(macroscopic)  equations  are  invalid  and  must  be  replaced  by 
Boltzmann's  (microscopic)  equation.  In  the  interest  of  simplicity, 
however,  we  shall  disregard  this  fact  and  assume  that  the  quasi- 
hydrodynamical  equations  are  always  valid.  We  must  remember 
though  that  the  predictions  derived  from  the  latter  equations  are 
only  qualitative. 

The  relevant  equations  are  Maxwell's  equations  and  the  equations 
of  motion  of  the  electrons,  ions,  and  neutrals.  21  (The  mks  system  of 
units  will  be  employed. ) 


V  xE=  -poH  =  (5) 


V  X  H=  J  +  =  ?  +  k^f 


=  jwP  +  jkowS  ,  (6) 


where 


T=  Ti+72=  NjCj?!  +  N,e2?,  =  +  P 


P=  Niejri  - 


(7) 

(8) 


"^1^1  ■  ®1 


[  -  rj2(rj  -7^) 


-  ^13^^  -  *•3’ 


(9) 


8 


[  gt  -  r2i(r2  -  Tj)  -  FjjCfj  -  Tj)  (10) 


•  • 


'"3^3=  -r3i(r3-r,).  r  32(73-^2'  . 


(11) 


where 


ri2=  +  m^)  .  etc. 


In  the  above  equations  and  throughout  this  report,  the  subscripts 
1,  2,  and  3  refer  to  electrons,  ions, and  neutrals,  respectively.  Thus, 
®1  ~  ■®'  ®2  “  ®  ®3  ~  where  e  =  1.  60  x  10*^^  coulomb.  The 

dot  means  partial  derivative  with  respect  to  time,  and  as  we  are 
assuming  d/dt  ^  3/8t  and  a  variation  with  time  exp  (jot),  the  dot 
means  multiplication  by  jc*).  In  writing  (9)  and  (10),  we  have  assumed 
that  the  omitted  terms  are  negligible  by  virtue  of  conditions  (e)  and  (f). 
Because  of  conditions  (g)  and  (h),  eqs.  (9)  and  (10)  can  be  simplified  to 

"•I'l  =  =1  [  ^  +  ■  rjjrj  (9a) 

and 


.  63  [  E  t  (j'2  I'o”’]  ■  '■  23"2  '  <“>*> 

It  should  be  noted  that  in  (9)  and  (10).  H  is  the  magnetic  field  of  the 
earth  alone,  for  the  magnetic  field  of  the  wave  is  negligible  compared 
with  it, 

Consider  a  right-handed  rectangular  coordinate  system  with 
the  x-axis  pointing  vertically  upward,  the  y-axis  perpendicular  to  the 
magnetic  meridian,  and  the  z-axis  in  the  magnetic  meridian  so  that 
the  earth's  magnetic  field  lies  in  the  xz  plane.  Simple,  but  rather 
laborious  algebraic  majiipulation_s  of  (9a),  (10a),  and  (7)  then  yield^B 
J=  ^  where 
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“1-4  j“l^U  -^Ix^lz 


Qj(q 


2 

1 


■  ^Ix^z  ‘■’“l^x  ■  ^z 

(12) 

and  a  similar  equation  for  J2  with  the  subscript  2  replacing  the 
subscript  1.  where 


=  Ne^/m.k 
0^  '  1  o 


Wj  =  e^B/m^ 


Cl  =  UUj/u 


Qj  =  u^d  -  jVj3/w)/«^ 


1 


Nj  -  =  N  =  electron  and  ion  number  density. 

(Quasi-neutrality  is  assumed.  ) 


(13) 

Consider  a  plane  MHD  wave  which  is  excited  in  the  ionosphere 
and  which  is  obliquely  incident  at  the  lower  edge  of  the  ionosphere- 
atmosphere  interface,  but  with  its  plane  of  incidence  in  the  magnetic 
meridian  or  the  xz  plane  as  shown  in  Fig.  2.  Call  0^  the  angle  df 
incidence  at  x  =  80  km  and  ^  the  angle  of  refraction.  From  Snell's 
Law, 


10 


sin  ^ 
sin  0, 


A, 


c 

A, 


=  n, 


(14) 


or 


sin  ^  =  n^  sin  0^  =  sin  0 


(15) 


In  (15)  n  is  the  value  of  the  index  of  refraction*  at  the  height  x  in  the 

ve"rti>^r^*  ®  makes  with^he 

vertical,  as  shown  in  Fig.  2. 

As  the  wave  propagates  downward,  all  field  quantities  will  vary 
as  explj^t  -  ---(sin  9a  -  cos  9..,)  i„  Region  II  and  as  exp[j(„t  -  k  ^ 
sin  dz-  cos  dx)l  in  Region  I,  where  k.,  =  2nf/c  is  the  vacuum  waJe- 

,;"®‘5in  d  V"®n‘cos’  9“P,V“^',-,  T  ^  ^  -  cos  9  x))  =  exp[j„ 

—  and  in  taking  partial  derivatives  of  the  field 

quantities  in  the  ionosphere,  we  may  set  V  =  0,  V  =  -  J^sin  j 

-  jk^sin  i  =  .  jk^cos  Z.  where  Z  is  the  angle  the  Refracted  ray  makes 

with  the  z  axis.  In  the  atmosphere,  we  have  V  =jk^cos<i.  V  =  0 

and  V  =  .jk  sin  <5  .  x  ^  v  .  y  u 

z  v 

meridian  InH  the  plane  of  the  magnetic 

that^''  respectively,  then  a  similar  argument  shows 


'^y  -  ^  •  ^.7  =  COS  Z 

/  Z  y 


in  Region  II 


and 

(16) 


V  =(jk^  cos  X.  -jk^  cos  Y.  -jk^  cos  Z)  in  Region  I. 


One  of  the  important  effects  of  the  inclusion  of  collisions  is  that  n  is 

complex  and  the  waves  are  attenuated  as  they  propagate.  In  general 
0  and  Ol  are  complex  angles.  f  t'  ^  .  m  general. 
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X 


REGION  m  EXOSPHERE 


X:550  KM 


REGION  H  IONOSPHERE 


\\\\\\\\\\\\\\\\\\\m 

PLANE,  INFINITELY  CONDUCTING  EARTH 


Fig.  2.  Plane  MHD  wave  excited  in  the  ionosphere  and  incident 
obliquely  at  the  lower  edge  of  the  ionosphere  and  being 
reflected  and  refracted.  The  plane  of  incidence  is  in 
the  magnetic  meridian.  {Not  drawn  to  scale.  ) 
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«motDh?r'.'  “<i  emerge,  in  the 

w.Tl  .  u  ■  “"?**  “■“*'  ef  i‘  '>y  the  uee  of 

wJve\  in7homo".'.n“.V  ‘h*  PropegMion  of  plane 

Chlnf.l.  e  medium.  A.  a  uaeful  reference  we  may  cite 

..c«ln  1  ^  u\°‘  “  'he  »et  of  ihi. 

.  j  concerned  mainly  with  etudying  the  hietory  of  the 

incident  and  reflected  MHD  wave  in  the  iono.phere 


By  ueing  (7),  (12),  (13),  and  (16)  in  Maxwell',  equation. 


and  (6).  we  obtain 


(5) 


-  jk^co.  Y  t  jk^co.  Z  Ey  = 


-jk^co.  Z  E^  -  V  -juPoHy 


+  jk  COS  Y  E  =  -jwu  H 
X  y  V  X  ■'^^0  z 


(17) 

(18) 
(19) 


and 


-  jk  cos  Y  H  +  jk  cos  Z  H 

V  z  ■'  V  y 


■lx 


“ifi  - 


E 

X 


/  ^Ix^lz  ^  ^2x  hz 

\“l  <  “l  -  ^1  )  ®2(  4-^1) 


E 

z 


(20) 


13 


■,\cos  Z 


-  (2) 


j 


( 


’Iz 


-  c 


(21) 


"y  *  Hx 


I  5,xC,x 


E 

X 


+ 


+  (2) 


E  + 
y 


“i  -^iz 

]  2 - ZT 

“l  <  ^ 


u  meant,  h  our  treatment  •'>’» 
looting  and  that  ts  why  ”  ,  ;  toc  i"’ 

K:  “kS  '•''••">  •■=;■;.“■'* I'. 

and  can  be  omitted  in  the  MHD  app^ox'imatiom®  displacement  current 
quantities,  we  obtain^he  desired  result. 


'y.  ^z.  riy.^  and  Doing  this  and  defining  the  ffllowine 

T*OPll1f  «r«p.  /  ^  \  ® 


viz.  (26). 


Q.  = 


1  z 


+  (2) 


Q.  = 


J^lx 
“2 — 
"i  -  ^1 


+  (2) 


Q3. 


1  - 


(2) 


Qa  = 


Ui  (  Uj 


(23) 


^Ix^lz 

r~ 


+  (2) 
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“l  • 

°5-‘— 7  i  2,  HZ)  . 


I  2  ,2i 

a,(  a  -C  ) 

~7~2 — — rz — Z~  • 

“l*  »1  -  Cl)  -  (“i  -CiJ 


°7  =  Q6°4  • 


°8  =  06°1 


“9  = 


)(«,  -  1)C,^ 


Vv?' 


+  (2)  .  Q9  *  Qg  +  Q 


1  ’2 


ja,Q 


'10 


1-9^  jQ2Q9 

rrr  "  -r 


0..4  =  1  + - ^  ^ 


lx 


2x 


Qll=  1 


^  "  •’^1x^9, 


1  -  JC2XQ9. 


(23) 


Qy  cos  Z 


M  =  jk 

V 


Qy  cos  Y 


Qg -  cos  Y  cos  Z 
Q|j  -  cos^  Y 


-Qg  cos  Z 
-Qg  cos  Y 
-Qjq  +  cos^  Z 
Qg  +  cos  Z  cos  Y 


Qg  cos  Z  cos  Y 
-1  +  Qg  cos^  Y 
Qg  cos  Y 
-Q^  cos  Y 


1  -  Qg  cos^  Z 
-Qg  cos  Z  cos  Y 
-Qg  cos  Z 
Qy  cos  Z 

(24) 
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E 

z 


f 


E 

y 


0  y 


where 


'’0= 


=  377  ohms. 


V  F  =  M  ■  f 

X 


(25) 


(26) 


,  Clemmow^.  4  terms  the  coupled  form  of  the 

differential  equations  governing  the  propagation  of  MHD  waves  in  the 
ionosphere.  If  we  are  to  trace  the  path  of  an  MHD  wave  excited  in 
the  mnosphere.  we  must  solve  (26)  subject  to  the  following  boundary 


atx=0 


(e^  X  E)  and  (e^  x  H)  continuous  at  x  =  80  km  (27) 


^  approaches  infinity,  we  have  only  an  outgoing 
plane  wave. 

In  addition  Table  I  shows  how  the  properties  of  the  ionosphere  vary 
with  altitude  or  x  (see  Fig,  2). 
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Except  for  very  special  cases  (such  as  (a)  vertical  incidencell 
when  cos  Y  =  cos  Z  =  0;  (b)  oblique  incidence  and  vertical  magnetic 

field  when  4,  ;  =  0;  and  (c)  propagation  in  the  magnetic  meridian 
and  at  the  equator  when  cos  Y  =  0,  cos  Z  =  1,  ;  =  0.  r  =  r),  it  is 
well  nigh  impossible  to  solve  (26).  In  our  problem,  we  have  the 
additional  complication  that  the  geometries  optics  approximation  so 
widely  used  at  radio  frequencies  or  microwave  frequencies  are  invalid 
because  the  properties  of  the  ionosphere  change  so  rapidly  in  a  wave¬ 
length.  The  last  statement  can  be  verified  by  remembering  that 
X  S?  100  km  and  by  studying  Table  I  and  Fig.  3. 

^  *‘®^crences  3  and  4.  Clemmow  gives  a  very  detailed  and  clear 
solve  (26)  in  the  three  cases  mentioned  above  with 
me  aid  of  a  modern  high-speed  computing  machine  such  as  the  IBM  709. 
The  general  method  will  be  indicated  here,  but  for  details  the  interested 
reader  is  referred  to  Clemmow's  papers. 

The  matrix  M.  eq.  (24).  is  a  function  of  x.  Presumably  it  is 
possible  in  any  small  region  of  the  ionosphere  to  find  a  local  coordinate 
system  in  which  many  of  the  elements  are  zero.  When  this  is  the  case 
^e  solution  of  (26)  is  greatly  facilitated.  Instead  of  the  column  vector ' 
t  of  (25),  imagine  that  the  field  quantities  are  express^  in  terms  of 
another  column  vector  G  related  to  f  by  the  matrix  0.  Thus. 


f  =  U  .  C  ,  (28) 


^  Assuming  that  U  is  nonsingular  so  that  it  has  an  inverse 

1  /  'u  substituting  (28)  into  (26)  and  multiplying  on  the 

left  by  U*l, 


-  (U-IMO)  .  G  =  -U-l  (V^U)  •  G 


If  L  =  (U-IMU)  is  a  diagonal  matrix,  the  solution  of  (29)  is  con¬ 
siderably  easier  th^  the  solution  of  (26).  From  a  well-known  theorem 
in  matrix  theory,  L  is  a  diagonal  matrix,  provided  the  roots  of  the 
characteristic  equation, 


det 


0 


(30) 


are  distinct.  Then 


17 


ALTITUDE  ,  KM 


Fig.  3.  Electron  and  neutral  particle  number  densities 
in  the  ionosphere  and  exosphere.  (From 
reference  Z7.  ) 
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The  equations  to  be  solved  are 


(31) 


-  L  .  ^  =  -U-1  (v^U)  •  G 


(32) 


where 


2  =  u-1  .  ? 


(33) 


mentioned  in  Section  I.  Francis  and  Karplus^^  have  solved 
(32)  for  vertical  incidence  with  the  boundary  conditions  prescribed  by  (27), 
the  geometry  given  by  Fig.  2.  and  the  properties  of  the  ionosphere  given 
by  Table  I.  Of  the  many  results  they  obtained  for  this  simple  case,  those 
that  interest  us  most  are  the  (complex)  transmission  and  reflection  coef- 
ficients  as  a  function  of  incident  angular  frequency  given  by  their  Table  I,  25 
Although  the  equations  contained  in  (32)  are  much  more  complicated  than 
those  studied  by  Francis  and  Karplus,  their  paper  and  Budden's2  are  worth 
careful  study  if  it  is  deemed  important  to  solve  (32). 

Now,  a  system  of  four  linear  differential  equations  is  equivalent 
to  a  single  fourth-order  equation  which  has  four  independent  solutions. 

Near  the  top  of  the  ionosphere,  two  of  these  are  upgoing  waves  and  two 
are  downgoing  waves.  24  There  are  therefore  two  independent  solutions 
which  satisfy  the  radiation  condition  of  (27).  One  workable  procedure  is 
t^o  start  with  one  "initial”  solution  and  integrate  (32)  by  means  of  a  step- 
by-step  process  downward  through  the  ionosphere  until  a  point  below 
the  ionosphere  is  reached.  The  resulting  solution  gives  the  total  values 
of  G  below  the  ionosphere.  It  is  then  separated  into  an  upgoing  and  a 
downgoing  wave,  each  being,  in  general,  elliptically  polarized.  The 
ratio  of  the  amplitudes  of  these  two  waves  gives  a  reflection  coefficient, 
but  this  applies  only  to  an  incident  wave  of  a  particular  elliptical  polariU- 
tion,  which  is  of  no  special  interest.  The  process  is  repeated  with  the 
other  "initial"  solution,  and.  in  general,  the  resulting  upgoing  and  down- 
going  waves  below  the  ionsphere  have  different  polarizations  from  the 
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first  pair.  One  then  takes  proper  linear  combinations  of  the  two  solutions 
such  that  the  incident  wave  of  the  combination  is  of  unit  amplitude  and 
plane  polarized  either  in  the  plane  of  incidence  or  perpendicular  to  it. 

The  associated  downgoing  wave  thengives  two  of  the  components  of  the 
reflection  coefficient.  The  exact  procedure  which  must  be  followed  to 
determine  the  reflection  coefficients  is  clearly  described  in  paragraph  8 
of  reference  2. 
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III.  THE  IMPORTANCE  OF  NEUTRALS  IN  DETERMINING  THE 

NATURE  OF  THE  MHD  WAVES 

From  Fig.  3  and  Table  I,  it  is  seen  that  throughout  the  entire 
mnosphere  the  number  density  of  the  neutrals  greatly  exceeds  the  number 
density  of  the  charged  particles.  An  immediate  consequence  of  Nn 
greatly  exceeding  Ng  and  Ni  is  that  in  Table  I.  v  and  v-  are  the 
electron  and  ion  collision  frequencies  with  the  neutrals  alone,  for 
'^ei'  Hei^ndv^  are  all  negligible  compared  with  v  and  v. 

en  in* 

A  very  basic  and  important  property  of  a  plasma  is  its 
tendency  toward  electrical  neutrality.  The  theory  predicts  and  measure- 
ments  confirm  that  quasi-neutrality  or  Ng  ^  N^  holds  in  a  plasma  in 
equilibrium.  A  good  measure  of  the  correctness  of  the  latter  approxima¬ 
tion  IS  the  smallness  of  the  Debye  length  ij)  compared  with  all  relevant 
macroscopic  dimensions  such  as  the  wavelength,  mean  free  path, 
cyclotron  radius,  and  the  dimensions  of  the  system.  Although  the 
temperature  as  a  function  of  altitude  is  not  an  accurately  known  property 
of  the  ionosphere,  the  temperature  can  be  estimated  within  reasonable 
limits.  It  is  instructive  to  note  that  the  maximum  value  of  the  Debye 
length  throughout  the  ionosphere  is  approximately  0.  5  cm. 

In  a  fully  ionized  plasma  the  velocity  of  HM  waves  is 


where  p  is  the  mass  density  of  the  plasma  p  =  pg  +  pj.  In  a 
partially  ionized  plasma,  the  question  naturally  arises:  "What  is  p 
in  eq.  (34)?"  From  their  definitions,  it  is  obvious  that  if  v;<<Cr  u-, 
the  motion  o_f  the  ions  is  constrained  by  the  magnetic  field  to  be  ^ 
parallel  to  B.  When  v  is  very  large  compared  with  the  collision 
frequencies  of  the  charged  species  with  themselves  and  with  each 
other,  the  neutrals  in  colliding  with  the  ions  effectively  drag  them  in 
the  direction  parallel  to  B.  Before  the  latter  statement  is  true, the  mass 
of  the  ions  and  the  mass  of  the  neutrals  must  be  approximately  equal 
so  that  considerable  momentum  can  be  exchanged  when  an  ion  and  a 
neutral  collide.  As  we  can  safely  assume  that  mj  =  m^,  it  follows 
that  in  the  lower  ionosphere 


B(B  .  Vn) 
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electron^'^^nlt^^a  (10a)  and  (11)  are  the  equations  of  motion  of  the 
electrons,  ions,  and  neutrals,  respectively.  Using 


we  can  write  (11)  as 


m  V 

n  ■"  m  +  m  ^'^n  '  '^e^ 
n  e 


m.v  . 

1  ni  -  - 

m  +  m.  ^'^n  "  '^i^ 
n  1 


The  exact  expressions  for  v  ariH  v  .  fUo.  r 

between  neutrals  and  electrons  and  n?Strals  and'ions,  respecdvefy '^are^^* 
?f  we  uTthreLo®  temperatures  and  densities  typical  of  the  ionsphlre. 

JuZn;,\:T:rV:rTZi': 


=  [  2wKT(m^+m^)/,„,_^„^,j  1/2.. 


N  v  /N 
e  en'  n 


w^'he  ‘'"I-  "collision  distance" 

or  the  distance  of  closest  approach.  ) 

V  =  2orv°"’l^^“'v‘'  “ 'oc  sli'io'ios  between  80  and  200  km. 

®o  V  '  '  '^ne  "  ^00  V  When  we  estimate  m„  =  2  x  lO'l  m 

and  remember  that  m„  ^  m^^the  first  term  on  the  right-hand  sideT' 

ind’  wltrifo^^hordLT  W^retcett;^”"''  ^PP-*‘-atinn 


^  -  1  +  2ja;/e_^.  =  i  t  2j„Ny(N.  v.)  e  ]  +  jQ  .  ,38, 

of  Q  ■"  (38)  relative  to  one,  we  can  draw 
the  following  conclusions:  (a)  Q>.1;  i„  ,h,s  case  the  neutrals  are 

fb?  q‘Ss^ Virjhr''' 90  degrees  out  of  phase. 

lith°htloUon  oJ  the  mns  i”'’.':™  “''"ohronired 
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V.  at  V 
1  n 


(39) 


When  (39)  holds,  we  must  set 
under  certain  conditions,  the 
determining  Vy^. 


P  -  Pe  Pi  Pn  (34).  In  other  words, 
presence  of  the  neutrals  is  important  in 


I 

I 
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glossary 


In  this  report  the  mks  system  of  units  is  used, 
important  symbols  are  listed  below. 


Some  of  the  more 


magnetic  induction 


speed  of  light  in  vacuo  »  3x10®  meters/sec 
speed  of  sound 

e -folding  distance  of  MHD  waves,  eq.  (2) 
electric  displacement 
electric  intensity.  5  = 

e  as  a  subscript  means  electrons,  magnitude  of  electron 
charge 

plasma  frequency  =  9  x  10^  cps 

frequency  of  the  MHD  or  EM  wave,  Zrrf  =  w 
total  magnetic  intensity 

static  magnetic  intensity,  earth's  magnetic  intensity 
perturbed  magnetic  intensity.  H  =  llTl/li?  |«^1 

unit  matrix,  eq.  (30) 
i  as  a  subscript  means  ions 
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/ 


j 

J 

k 

k 

( 

K 

L 


m 

e 


mi 


N 

e 

N. 

1 

N 

n 


n 


n 


P 


the  imaginary  unit.  j2  =  .1.  Field  quantities  vary  with 
time  as  exp  (jwt). 

current  density 

X  =  wavelength  of  MHD  waves 

permittivity  of  free  space  =  8.  854  x  10-12  farad/meter 

Boltzmann's  constant 

characteristic  length";  LgQ  is  the  characteristic 
length  at  an  altitude  of  80  km. 

Debye  length  =  6.  .90(T/n  )  cm 

e 

electron  mass,  also  designated  by  mj 
ion  mass,  also  designated  by 

electron  number  density  per  cm^.  also  designated  by  Nj 

ion  number  density,  also  designated  by  N2 

neutral  number  density,  also  designated  by 

unit  normal  at  the  boundary  of  a  discontinuity,  n  as  a 
subscript  means  neutrals 

index  of  refraction  =  c/V  . 

ph 

pressure 
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polarization  vector,  eq.  (8) 
position  vector,  (x,  y,  z) 
temperature  in  degrees  Kelvin 
Alfven  velocity,  eq.  (4).  Fig.  1 
phase  velocity  =  u/k 

mean  fluid  velocity  in  an  element  of  volume  AV 
mean  electron  velocity  in  an  element  of  volume  AV 
mean  ion  velocity  in  an  element  of  volume  AV; 


P;  +  P 
i  '^e 

average  ion  molecular  weight,  Table  I 

X  coordinate,  the  altitude 

■  j 13/w)/wqj  .  eq.  (13) 

ratio  of  the  sum  of  the  gas  and  plasma  pressure  to  the 
magnetic  pressure,  eq.  (3) 

angular  frequency  of  the  MHD  or  EM  wave,  u  =  Zirf,  eq.  (7) 

electron  cyclotron  frequency  =  eB/mg  =  1.  8  x  lO^B 
(rad/sec,  gauss) 
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U).  = 


w. 


ion  cyclotron  frequency  =  eB/m- 

pla.ma  frequency  =  V4wN^eVm^  =  5.  5  x  10^  VnI  rad/sec 


wavelength  of  MHD  waves 

magnetic  permeability  =  1.  26  x  lO’t  henry/meter  B  =  u  H 

^0 

mass  density 


v . 


electron  collision  frequency,  v  =  v  +  w  .  v 

e  ee  ei  en 

ion  collision  frequency 
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1-type  particles  and  Z-type 
particles,  is  similarly  defined 


12  ”^1*^2  ^1  ^2^ 
scalar  conductivity 


tensor  conductivity,  eq,  (12) 


uu./uo 


V,,.  Vj 


y  z' 


1.  2.  3 


neutmTs  ‘,V  ^  •  mns,  and 


27 


REFERENCES 


1.  Gershman,  Ginzburg  and  Denisov,  Uspekhi  Fiz.  Nauk.  6l  , 
561-612  (1957);  AEC-tr.3493. 

2.  K.  G.  Budden.  Proc.  Roy.  Soc.  A  227,  516-537  (1955). 

3.  P.  C.  Clemmow  and  J.  Heading,  Proc.  Cambridge  Phil.  Soc. 

Pt.  2,  319-333  (1954). 

4.  K.  G.  Budden  and  P.  C.  Clemmow,  Proc,  Cambridge  Phil,  Soc. 
53,  Pt.  3,  669-682  (1957). 

5.  P.  S.  Epstein.  Nat.  Acad.  Sci.  _1^.  627-637  (1930). 

6.  L.  M.  Brekhovskikh,  Waves  in  Layered  Media  (Academic  Press, 
Inc.,  New  York,  19607i 

7.  S,  Akasofu.  J.  Atom.  Terr.  Phys.  _1_5*  156-160  (1959), 

8.  B.  Lehnert.  Tellus  8.  241-251(1956). 

9.  S.  Akasofu,  J.  Atom.  Terr.  Phys.  _1_8,  160-173  (I960), 

10.  A.  J.  Dessler.  J.  Geophys.  Research  M,  397-401(1959). 

11.  W.  E.  Francis  and  R.  Karplus,  J.  Geophys.  Research  65, 
3593-3600  (1960). 

12.  S.  Matsushita.  J.  Geophys.  Research  M,  1149-1161(1959). 

13.  R.  Liist.  Fortschritte  der  Physik  7.  503-558  (1959);  see 
especially  p.  528. 

14.  H.  Alfven,  Cosmical  Electrodynamics  (Oxford  University  Press, 
London,  1950);  see  especially pp.  78-88. 

15.  G.  H.  A.  Cole.  Advances  in  Physics  452-497  ( 1956);  see 

especially  p.  462.  ~ 

16.  J.  A.  Fejer,  J.  Atom.  Terr.  Phys,  _18.  135-146(1960). 

17.  J.  W.  Dungey,  Penn,  State  Univ,  Scientific  Report  No.  57,  1954. 

18.  J.A.  Fejer,  op.  cit.  ,  pp.  137-138. 


19.  L.  Spitzer,  Jr.  ,  "The  Terrestrial  Atmosphere  above  300  km,  " 
^he  Atmosphere  of  the  Earth  and  the  Planets  (University  of 
Chicago  Press,  Chicago,  Illinois,  1952),  edited  by  G.  Kuiper. 


28 


20. 

21. 

22. 

23. 

24. 

25. 

26. 

27. 

28. 


W.  Viezee,  Paper  1972  of  the  RAND  Corporation,  April  18,  I960. 

A.  Schluter.  Zeitschrift  fiir  Naturforschung  i^,  73-78(1951). 

J.  Stratton,  Electromagnetic  Theory  (McGraw-Hill  Book  Company. 
Inc..  New  York.  1941).  see  especially  pp.  327-330. 

J.  Ratcliffe,  The  Magneto-Ionic  Theory  (Cambridge  University 
Press,  LondorT  1959),  see  especially  Chap,  18. 

H.  G.  Booker,  Proc.  Roy.  Soc.  A  155.  235-257  (1936). 

W.  E.  Francis  and  R.  Karplus,  op.  cit.  .  p.  3595. 

S.  Chapman  and  T.C.  Cowling.  The  Mathematical  Theory  of 
Gases  (Cambridge  University  Press,  London, 

1 tSo ) , 

F.S.  Johnson.  J.  Geophys.  Research^,  577-584(1960). 

Dessler,  Francis, and  Parker,  J,  Geophys.  Research  65 
2715-2719  (1960).  — ’ 


29 


/ 


UNCUSSIFIED 


UNCLASSIFIED 


